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Abstract. We present a new time-dependent inhomogeneous jet model of non-thermal blazar emission. 
Ultra-relativistic leptons are injected at the base of a jet and propagate along it. We assume continuous reac- 
celeration and cooling, producing a relativistic quasi-maxwellian (or "pile- up") particle energy distribution. 
The synchrotron and Synchrotron-Self Compton jet emissivity are computed at each altitude. Klein-Nishina 
effects as well as intrinsic gamma-gamma absorption are included in the computation. Due to the pair 
production optical depth, considerable particle density enhancement can occur, particularly during flaring 
states. Time-dependent jet emission can be computed by varying the particle injection, but due to the sensi- 
tivity of pair production process, only small variations of the injected density are required during the flares. 
The stratification of the jet emission, together with a pile-up distribution, allows significantly lower bulk 
Lorentz factors, compared to one-zone models. Applying this model to the case of PKS 2f 55-304 and its 
big TeV flare observed in 2006, we can reproduce simultaneously the average broad band spectrum of this 
source from radio to TeV, as well as TeV light curve of the flare with bulk Lorentz factor lower than f5. 

1 Introduction 

It is widely admitted that the blazar phenomenon is due to relativistic Dopplcr boosting of the non-thermal 
jet emission taking place in radio-loud Active Galactic Nuclei (AGN) whose jet axis is closely aligned with the 
observer's line of sight. Blazars exhibit very broad spectral energy distributions (SED) ranging from the radio to 
the gamma-ray band and dominated by two broad band components. In the Synchrotron Self Compton scenario 
(SSC), the lowest energy hump is attributed to the synchrotron emission of relativistic leptonic particles, and 
the highest one is attributed to the Inverse Compton process (IC) of the same leptons on the synchrotron photon 
field. Broad band observations of these objects are crucial to understand the jet physics and to put reliable 
constraints on jet parameters. 

The most extreme class of blazars are the highly peaked BL lac sources (HBL), where the synchrotron/Inverse 
Compton components peak in the UV/X-ray/gamma-ray (GeV up to TeV) range. These objects are well known 
to be highly variable in all energy bands. Perhaps the most extreme example of this extraordinary variability 
behaviour has been caught by the HESS instrument with the big flare of PKS 2155-304 during summer 2006 
(Aharonian et al. 2007). 

In this case, the observed variability time scale {'^ 200 sec) in the TeV range implies a minimum bulk Lorentz 
factor greater than 50 (Begelman, Fabian, & Rees 2008) assuming an homogeneous one zone model. However, 
such high values of the bulk Lorentz factor are in contradiction with constrains derived from other observational 
evidence (Urry & Padovani 1995; Henri & Sauge 2006 and references therein). Furthermore, one-zone models 
are unable to fit the entire spectrum, the low energy radio points being generally attributed to more distant 
emitting regions. We present here a new approach, unifying small and large scales emission regions: we consider 
that the radio jet is actually filled by the same particles originating from the high energy emitting region, at 
the bottom of the jet, that have propagated along it. We describe thus the emitting plasma by a continuous 
(although variable) particle injection, submitted to continuous reacceleration and radiative cooling. This model 
fits well into the two-fiow framework originally proposed by Pelletier (1985) and Sol et al. (1989) where a non 
relativistic, but powerful MHD jet launched by the accretion disk, surrounds a highly relativistic plasma of 
electron-positron pairs propagating along its axis. The MHD jet plays the role of a coUimater and an energy 
reservoir for the pair plasma, which is responsible for the observed broad band emission. 
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Fig. 1. Sketch of the jet geometry. See text for the 
signification of the different parameters. 



2 Description of the model 

2.1 Geometry of the model and particle energy distribution 

We consider that the relativistic plasma propagates in a stationary funnel (associated to the last magnetic 
surface of the surrounding MHD jet) that we parametrize as a shifted paraboloid shape (power-law with an 
index oj). The magnetic field scales with the jet radius to the power — A. We consider that the jet is continuously 
accelerating, starting from rest {Tb = 1 at z = 0) and reaching an asymptotic value T^oo on a typical scale Zq. 

The MHD jet that surrounds the radiative pair plasma transfers its energy to the pairs via second order 
Fermi acceleration mechanism (particle- wave interaction) . As a result of this acceleration mechanism, the energy 
distribution function (EDF) of the electron-positron plasma is assumed to be a relativistic maxwellian (or "pile- 
up") distribution (Sauge & Henri 2004 and references therein). As they propagate inside the structure, the 
particles lose energy via synchrotron and inverse Compton cooling. Particles are continuously re-heated by the 
surrounding MHD jet, via a Fermi H acceleration process. The acceleration rate Qacc{z) is parametrized as a 
power-law with an index — C and a normalisation factor Qq. To avoid energy divergence, the heating is stopped 
after an altitude Z^- The relativistic maxwellian's "temperature" 70 is determined by balancing radiative cooling 
and re-heating. Because of absorption of 7-ray photons via the pair creation process, new particles are created 
inside the jet, which yields to the increase of the particle flux '^{z,t) = J n{'-f, z,t)S{z)rtf3bcd'y, where S{z) is 
the surface of radial section of the jet. This quantity is computed via a continuity equation that takes into 
account the pair creation term. Fig. 1. displays the different parameters used in the model. More details about 
the parametrisation we use can be found in Boutelier et al. (2008). 



2.2 The jet spectrum 

We compute the emissivity at each altitude in the jet by assuming a SSC process for the radiative mechanism. 
The total intensity of the jet is then determined by integrating the emissivity all along the jet. Once injected at 
the base of the jet, the particles contribute first to the high energy part of the jet SED. As they propagate, their 
emissivity peaks progressively at lower energy, producing the low energy part of the spectrum. The spectral 
shape of the whole SED is not controlled by the local particle energy distribution (which is always locally a 
narrow pile-up), but rather by the z-dependencies of the jet radius, the magnetic field, and the acceleration rate. 
A constant injection rate would lead to a stationary emission pattern, which would be rather easy to fit. In 
reality, the observed instantaneous spectra are a complicated convolution of the whole history of the jet, keeping 
the memory of the whole past injection pattern. Boutelier et al. (2008) describes a procedure to extract the 
physical parameters of the jet from observed spectra, despite the fact that they do not correspond to a simple 
steady-state of the jet. 
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Table 1. Model parameters of the flaring and quiescent state. During the flaring state, the flux of injected particles 
varies between the indicated minimum and maximum values, following the injection pattern displayed in Fig. 2. The 
other parameters remain fixed. Ri, i?o, Zq, Zc are in unit of 10^*cm and are displayed on Fig. 1. 
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Fig. 2. Left : Fit of PKS 2155-304 data. Filled dots: average archival data (see text). Empty triangles: average 
HESS data from the big flaring night. Empty diamonds: "fake flaring" state low energy points Shaded area: 
enveloppe of archival X-ray data from BeppoSAX, SWIFT and XMM-Newton. Dot-dashed line: best flt of the " 
average fake flaring" spectrum. Dashed line: fit of the quiescent (^average) spectrum. Solid line: example of an 
instantaneous simulated spectrum. Right: Upper panel: HESS light curve above 200 GeV superimposed with the 
model (solid line). Middle panel: time dependent particle injection function used in the simulation. Lower panel: 
predicted light curves in X-ray (dashed line, left y-scale) and optical (dot-dashed line, right y-scale). The dotted lines 
mark the maximum of the different bursts of the injection function. 



3 Application to PKS 2155-304 

We have tested our model to the big flare event observed by HESS in July 2006 in PKS 2155-304 . We 
apply the method described in Boutelier et al. (2008) to the SED of PKS 2155-304. The average spectrum of 
PKS 2155-304 in the low frequencies range (< 10^^ Hz), is constructed using archival data from the HEASARC 
archive website (http:/ /heasarc. gsfc.nasa.gov/docs/archive. html). Based on the activity detected by HESS, we 
estimate that the duty-cycle is around 10%. Then to construct the "fake flaring" spectrum of this source, i.e. the 
spectrum expected if the jet was always in a flaring state, we combine the HESS data with the radio-to-optical 
ones corrected by a factor 10 (dot-dashed line in left plot of Fig. 2). 

An average "fake flaring" spectrum is shown in the left side of Fig. [2] in dot-dashed line. The corresponding 
best fit model parameters are reported in Tab. [1] Interestingly we only need a Doppler factor Ff, = 15 which is 
significantly below the values of ~ 50 inferred from one-zone model (Begelman, Fabian, & Rees 2008). The best 
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fit parameters of the TeV quiescent state (derived from averaged H.E.S.S spectrum in Aharonian et al. 2005) 
are also reported in Tab. [1] and the solution has been overplotted in the left side of Fig. [2] (dashed Hne). 

To reproduce the observed H.E.S.S. light curve during the big-flare event, we use an injection function 
$(z = 0, t) that is the sum of five "generalized Gaussian" shape, similar to the analysis made in Aharonian et 
al. (2007). Before the flare, the flux of particles <i>(z = Q,t) is assumed to be a crenel function that oscillates 
between quiescent and flaring states in agreement with the source duty cycle. At the right top of Fig. [2l we 
have reported the HESS light curve and the simulated one (solid line). The agreement between the simulated 
and the reported HESS light curve is very good. The instantaneous spectrum overplotted in solide lines on the 
left side of Fig. [2] agrees with the radio to TeV spectrum observed during the burst. 

4 Discussion 

Our time-dependent inhomogeneous jet model succeeds in reproducing simultaneously the broad band (from 
radio to TeV) spectrum of PKS 2155-304 as well as the TeV light curve during the big flare event of July 
2006. The key idea of the method is to decompose the blazar spectrum in "quiescent", low luminosity states, 
and "flaring", high luminosity states. The high energy part of the spectrum, coming from small-scale inner 
regions, is assumed to be, at any time, in one of these pure states. On the other hand the low energy part is 
a convolution over a large scale of the past history of the jet : it is thus rather a time-averaged spectral state 
mixing quiescent and flaring states in proportions given by the source " duty-cycle" . Moreover we do not require 
two different populations of emitting particles like in other models (e.g. Katarzyhski et al. 2003) but simply a 
continuous (although variable) injection of particles at the base of the jet that propagate along the jet structure. 
Pair production plays an important role to amplify the initial variation, as can be seen with the variation of 
the particle flux along the jet during the flaring state (see Tab. 1 and the right side of Fig. 2 (middle)): the 
pair current is amplified by a factor 30 at the end of the jet, while the initial current varies only by a factor 2. 

The model can also predict light curves at different wavelengths. As an example, the x-ray (2-10 keV) and 
optical (V band) light curves expected during the TeV flare have been plotted at the right bottom of Figl2l The 
x-ray luminosity exhibits almost simultaneous variations but with a lower amplitude (about 5 times smaller). 
On the other hand, the optical light curve shows a very different behavior, increasing all along the flare. This is 
due to to the large size of the optical emitting region that plays the role of a low pass filter. Consequently, the 
optical luminosity integrates the recent past history of the jet. These results are compatible with simultaneous 
multiwavelength observations made during the "Chandra night" (Costamante, private communication). 

In the present work, we do not specify the origin of the variability. Obviously, other plasma parameters 
can vary in addition to the injection density Nq. Most likely, variability can be triggered by a change in the 
acceleration rate described by Qacc- In a- plausible scenario, long term (year scale) variability implying the 
succession of quiescent and active states could be attributed to variations in the accretion rate, wheras the 
short (minute-scale) flares would be attributed to the instability to pair creation that develops only when the 
initial particle density is close to a critical threshold. 
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